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Abstract 

Previously the quark and meson properties in a many quark system at h- 
nite temperature have been studied within effective QCD approaches in the 
Hartree approximation. In the present paper we consider the influence of 
the mesonic correlations on the quark self-energy and on the quark propaga- 
tor within a systematic l/N^- expansion. Using a general separable ansatz 
for the nonlocal interaction, we derive a selfconsistent equation for the 1/Nc 
correction to the quark propagator. For a separable model with cut-off form- 
factor, we obtain a decrease of the condensate of the order of 20% at zero 
temperature. A lowering the critical temperature for the onset of the chiral 
restoration transition due to the inclusion of mesonic correlations is obtained 
what seems to be closer to the results from lattice calculations. 
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I. INTRODUCTION 



QCD motivated effective theories are the most promising approaches to the low energy 
behaviour of QCD and the meson physics in terms of quark and gluon degrees of freedom and 
symmetries. Starting from chiral quark model Lagrangians a perturbative approach to the 
occurence of a chiral condensate below a critical temperature in mean field approximation 
is usually considered. Simultaneously, the pseudoscalar Goldstone boson, the pion, occurs. 
Perturbation theory can be formulated in l/Nc-, where Nc is the number of colors [|T|. The 
leading order is the Hartree approximation, results are reported in Refs. A more 

general approach, where a nonlocal instantaneous interaction is applied, has been presented 
in Refs. [3^TT||. A still open but very important question is the influence of mesonic degrees 
of freedom which are neglected in the Hartree approximation. These degrees of freedom are 
supposed to be dominant in the low temperature limit. For the NJL model, an effective 



1/Nc- expansion which accounts for the mesonic fluctuations has been considered in [|T2 



However, the set of diagrams for the self-energy in next to leading order considered in this 
reference was not complete. This has been observed in Ref. where also the role of the 
scalar iso- vector mesons in the 1/Nc - approximation was discussed at zero temperature. It 
was shown that in the l/N^- expansion the Schwinger Dyson equation for the quark self- 
energy is different from the gap equation for the quark condensate and has to be solved 



separately. A complete collection of diagrams in 1/Nc was given in Ref. [|14| and recently 
studied in the chiral limit niQ = at T = by Ref. At T = 0, effects of the order of 
10% — 20% have been obtained in these approaches, showing that mesonic fluctuations play 
an important role. 

In this work, we consider the influence of mesonic correlations on the quark condensate 
at finite temperature. It is expected that such a calculation beyond the Hartree level of 
description will lead to corrections to the temperature behaviour of the quark condensate 
since the medium allows for mesonic degrees of freedom. The relation of a generalized gap 
equation to the thermodynamical potential of a quark meson plasma has been considered 
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in Ref. [0. The present paper is a first step for a consistent description of a meson gas at 



finite temperature within a chiral quark model. 

The paper is organized as follows: In Section |I| the nonlocal chiral quark model is briefiy 
introduced, which is used in Section |T| to derive a generalized formula for the quark con- 
densate in 0{1/N^ expansion. In Section ^ we include dynamical fiuctuations into the 
self-energy and treat the scalar und pseudoscalar contributions within the pole approxima- 
tion. The numerical results for a calculation within the NJL model at finite temperature 
are discussed in Section 0. 

II. THE MODEL 

Our starting point is the chiral symmetric effective Lagrangian in the quark sector of the 
following general form 

^ = gi(p)(7/.P^ - mo)qi{p) + £int , (1) 
where the interaction term 

Ant = ~qi{Pl)^l2Q2{P2)K{pi,P2,Pl',P2')q2'{P2')^-V2'qi'{Pl') (2) 

is given as a nonlocal generalization of the current-current type interaction. Here the 
matrices Afg denote the decomposition into the color (c), fiavor (f) and Dirac (D) chan- 
nels. In this work we restrict us to scalar and pseudoscalar channels. Therefore we choose 
Ai2 = [Ic ■ 1/ • 1d]i2 and A^g = [Ic ■ tj ■ 275]i2- 

The gluonic degrees of freedom do not occur explicitely in this effective approach to the 
low-energy sector of QCD. They are assumed to form a condensate which is responsible for 
the nonperturbative character of the quark- quark interaction in this domain. We make the 
phenomenological ansatz of an instantaneous interaction kernel, which can be formulated in 



a covariant way |Tl[]. We employ here a separable form for the nonlocal 4-point interaction 
of the form 
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K{pi,P2,Pi,P2) = -j;f9{ — ^ — )9{ — ^ — )V-P2,p;-p^ • (3) 

The Nc- dependence arises from the Fierz transformation of the quark current-current in- 
teraction in the colour singlet channel considered here, see e.g. For our numerical 
calculations in Section |V| we use = 3. The dependence of the formfactor on the modulus 
of the three-momentum (|p| = p) has been discussed for different shapes, e.g. a Gaussian one 
{g{p) = exp[— (p/Acauss)^]) or the well-known NJL type interaction (^'(p) = 0(1 — p/Anjl)), 
see [0 . Note that the potential does not depend on the energy and we obtain therefore the 
NJL- model with a three-momentum cut-off. The spectral properties of the quark model 
defined by the Lagrangian (0) are obtained from the single particle propagator 

^12(^1^2) = [G(pi)lcl/]i2'5pi,p2 , (4) 

which is a diagonal matrix in color, flavor and momentum space. The matrix element G{p) 
obeys the Dyson equation 

G{p) = [G^\p)-j:{p)]-' , (5) 

where G'o^(p) = j^p^ — mo is the vacuum Green function, see Fig. 1. The self-energy 
S(j9) is defined by an analysis of all one particle irreducible diagrams contributing to the 
propagator. Having the single particle propagator at our disposal, the physical quantity 
of interest which is straightforwardly evaluated is the quark condensate. For our separable 
potential we introduce the nonlocal quark condensate as 

< qq >= N^Nf ^ ^(p)Tr [Gip]] , (6) 
p 

where Tr stands for the trace over the Dirac space only. The finite temperature investigations 
are performed using the Matsubara technique p8|-p0|, where po = ioJn with the fermionic 
Matsubara frequencies Un = {2n + 1)ttT and J2p stands short for TJ^nl dp/{2iT)^. In 
order to obtain estimates for the quark condensate one has to make approximations for the 
self-energy. 
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The first step towards a systematic investigation of tlie Dyson equation (H) is tlie self- 
consistent Hartree approximation, see Fig. 2, 

S^[p;G^] = -KoNfg{p)Y^gik)Tr[G''{k)] , (7) 

k 

wliich defines upon insertion in tlie propagator in Hartree approximation 

G^(A;) = [Go'(A;) -S^[k;G^]]"' . (8) 

The Hartree self-energy (0) is a Dirac scalar and appears as a mass term in the propagator, 

m-^(p) = mo - ^(p)^ < qq (9) 

with a momentum dependence due to the nonlocality of the interaction kernel (|^). The 
quark condensate in Hartree approximation is, cf. Eqs. (|^) - (0), 

<gg>^ =N,NjY.9{^)T^\G''{k) 

k '- 

for details see e.g. Refs. P-p|pi?|. In this approximation, the magnitude as well as the tem- 
perature dependence of the dynamical mass generation is determined from the condensate 
only. Note that the restoration of the chiral symmetry at temperatures above the critical 
one (Tc ~ 200 MeV, [jlO|) is governed by the Fermi distribution function of quarks in the 
medium, f{E{k)) = {exp [{E(k)) /T] + 1}^^, where the quasiparticle dispersion relation 



E{k) = ^k2 + (m^(k))2 (11) 

contains the momentum dependent Hartree mass (P). 

It is, however, questionable whether the Hartree approximation is appropriate for the 
description of the nonperturbative low energy region of QCD where free quarks should be 
absent due to confinement. Since mesonic correlations are supposed to dominate the low 
energy excitation spectrum of the quark matter system, one has to study their influence on 
the results obtained within the Hartree approximation. A systematic perturbation theory 
for strong interactions is however lacking. Instead, one resorts to an expansion of diagrams 
to orders which we will investigate in this work at finite temperatures. 



III. l/Nc- EXPANSION 



The self-energy in the selfconsistent Hartree approximation appears of the order 0[1] 
as the leading term in the l/N^- expansion, as can be seen from Eq. (|^. In order to 
improve this approximation, we will study next to leading order diagrams, i.e. 0[1/Nc]- 
contributions. Therefore we make the following ansatze for the self-energy and for the quark 
propagator 

= S^[p; G] + G] + 0[1/N'^] , (12) 

Gip) = G^ip) + ^JG{p) + 0(l/iVf ) , (13) 

where the corrections to the self-energy depend in the general case on the full Green function 
G{p) and on the 4-momentum p. The corrections to the self-energy ^^[p;^] are not yet 



specified and will be discussed in the following section. Using the l/N^.- approximation (|T2 
for T,[p; G] in the Dyson equation (|^), the l/Nc- expansion to the propagator is given as 



G{p) = {G^\p)-np;G] 



G,\p) - S^[p; G^] - ^ [S^[p; 5G] + G^] + 0{1/nS 



(14) 



Expanding the l/Nc contribution in the denominator and comparing with ([TsD, we obtain 
a selfconsistent equation for SG{p) in the form 

SGip) = G^(p)S^[p; 5G]G''{p) + G''{p)5E[p; G^]G^(p) . (15) 

Note, that this consistent l/Nc- expansion for the quark propagator is a new result of this 
paper. In particular, the first term on the r.h.s. of ([ISD has not been considered in some of 



the previous approaches, see ||12| , pT| . In order to get a closed expression, we use the fact that 
the functional dependence of the Hartree selfenergy on the l/Nc- corrections to the quark 
propagator 6G is known from Eq. (|^). After insertion of S^[p; 6G] on the r.h.s. of Eq. (p!5[), 
we obtain 
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J2g{p)Tr[SG{p)] = -K,NfJ2 9'{p)Tr 



j:g{k)Tr[5Gik)] 

k 



(16) 



l-J-(O) p 

where the scalar quark loop integral J°^(0) is defined in Appendix ^. The l/N^- expansion 
of the quark condensate corresponding to that of the propagator ([T3|) and the definiton of 
the quark condensate (Bf) reads 



< qq > = < qq +6 < qq > +0[1/N'^ 



The 1/Nc correction to the condensate is obtained in closed form using the result (|T6| 



6<qq>=Z-NfY. ^7(p)Tr 



(17) 



(18) 



with a prefactor Z = 1/(1 — J'^(O)) as derived in ( p!6D coming from the l/N^- contributions 
to the Hartree self-energy S^[p, see in Fig. 3. This prefactor leads to a considerable 
rescaling {Z ~ 4) which has first been pointed out in Refs. 1]T^|IS| for the NJL model at 
zero temperature and is obtained here for the more general case of a nonlocal separable 
interaction at finite temperature. 



IV. MESONIC CORRELATIONS 

Within the chiral quark model as defined in Section ||, the complete set of diagrams 
contributing in (9[1/A^c] to the self-energy is given in Fig. 4. The double line corresponds 
to the RPA - type partial resummation of the chain of bubble diagrams, where the quark - 
antiquark loop in Hartree approximation defines the polarization functions J'^{p — k) in the 
scalar and pseudoscalar channel {(f) = cr, vr), see Appendix 
The (9[1/A^c] self-energy contribution is given by 

|p + k| 



k 



l-J^'ip-k) ' ^ 'l-J^{p-k) 



(19) 



The denominators 1 — J'^{p — k) occur due to the resummation and thus strong correlations 
can be described. Note that the l/N^ self-energy is a dynamical quantity and has not 
yet been solved in its complexity. The most dramatic effect is the occurence of collective 
excitations in the quark - antiquark channel when ReJ'^(P = M^) = 1 (and ImJ'^(P = 
M^) = 0) which correspond to mesonic bound states. In what follows we restrict us to the 
consideration of bound states only and use an expansion of the polarization function at the 
mesonic poles (pole approximation) which leads to the introduction of meson propagators 
and meson- quark-antiquark form factors g^j^qq (see Appendix ^ 

1 - J'^(P) ^ Mj-P^ NfKo ■ ^ ^ 

The full treatment of the RPA approximation which contains bound and scattering states 
is possible for the separable interaction and will be regarded in an additional work. 

Using the expression (^, for the self-energy and the short notation with = 1^ 
and = i75, we obtain 



6 < qq > 



pqq 



1 - J'^(O) J (27r)3 
dpo f dko 



5(P) / 77^9 



(27r)2 



1 



2^ M^-ik-p)^ 
After evaluation of the Dirac trace 



Tr 



{p)T'f'G^ {k)T'f'G" {p) 



(21) 



1 



Mj-{k-py 

-4 



Tr 



G^(p)r*G^(A;)r^G^(p) 
m(p) ± m(k) 



(p2 _m2(p))(A;2 -m2(k))[(A;-p)2 -M2] 
M| - (m(p) ±m(k))2 



m(p) 



1 



(p2_m2(p))2[(A;-p)2-M2 

1 



(22) 



(p2 - m2(p))2(P _ m2(k))[(A; - pY - M2] (p2 _ ^2(p))2(p _ ^2(1^)) 

we perform the Matsubara summation (see Appendix ^) and obtain as the final result for 
the 1/Nc mesonic contributions (^Tj) to the quark condensate 



5 <qq>^ 



pqq 



9ip) 



d^k 



|p + k| 



1-J'^(0)y (27r)3^'^'7 (27r)3^ V 2 

2m^(p) f l - 2f{E{p)) f{E{p))[l - f{E{p))] 
E^p) { 2E{p) T 
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^ l-2f{E{k)) l + 2n{E4k-p)) 
2E{k) 2E^{k - p) 

[1 - f{E{p)) - f{E{k))] [1 + n{E{p) + E{k)) + njE^jk - p))] 
EHp)E{k)E4k - p)iE^ik - p) + Eip) + Eik)) 

X <j E\-p)[m^{p) ± m^{k)] + m^{p)[Ml - [m^(p) ± m"{k)]^] (23) 



E^{k-p) + 2E{p) + E{k) E{p)f{E{p))[l-f{E{p))] 
E^{k-p)+E{p) + E{k) ^ T 



+[m - -m] 



+ [i?(p) - -i?(p)]) } 



with the energies E(f,(k — p) = -^/(k — p)2 + M| and the bosonic distribution function 
n{E) = [exp{E/T) The upper sign holds for the scalar, the lower one for the 

pseudoscalar meson, respectively. The (9[1/Ai"c] contribution (^3]) consists of two parts, and 
the numerical analysis shows that the contribution due to mesonic correlations is dominated 
by the first one, i.e. 

^, 9^ f d'P . . f ,( \p + k\ \ m"{p)[l-2f{E{p))] 

o^qq^ ^ 1- J-(0)y (27r)3^^^V (27r)3^ V 2 J E'^{p) 
■l-2f{E{k)) l + 2n{E4k-p)y 



2E{k) 2E4k - p) ; ' ^^^^ 

which has a simpler structure than (|23|) . 

In order to compare our results with previous works we will now discuss the T = case. 
At zero temperature the Fermi and Bose distribution functions vanish. In the T = limit 
of Eqs. (|TOD and (^) we obtain 

. - 6 alqq f d^p , . [■ d?k ^(\p + k\\(m^{p) (I 1 \ 

[Ml - [m^(p) ± m"{k)f][E^{k - p) + 2E{p) + E{k)]m" {p) 



E^{p)E{k)E^{k - p){E^{k - p) + E{p) + E{k)Y 
[m^(p)±m^(k)] \ 

E{p)E{k)E^{k-p){E^{k-p)+E{p) + E{k))) ■ ^ > 

In the following section we present the numerical evaluation and discussion of the above 
1/Nc- corrections to the quark condensate. 
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V. NUMERICAL RESULTS AND DISCUSSION 



In section |T| we have introduced a general nonlocal interaction kernel in separable form. 
In order to compare the numerical results with previous approaches within the NJL model, 
we will restrict us in this paper to the discussion of a cut-off formfactor 

|p + k|\ ^/ |p + k|\ , , 

'^=0 1-^W • (26) 



2Anjl 

The chiral quark model with soft formfactors as, e.g. a Gaussian one, has been discussed in 
Refs. 

After fixing the parameters of the model as described in Appendix 0, we obtain for the 
quark condensate in the Hartree approxim ation < uu >^= —(250 MeV)^ and for the quark 
mass = 300 MeV in agreement with the well-known data of the literature 

In the next step, discussed in Section |^ we have included dynamical selfenergy contri- 
butions due to mesonic correlations. Compared with the Hartree term (p!OD where we have to 
solve a one-loop integral, the next to leading order contributions are two-loop integrals which 
after summation over both Matsubara frequencies {ko,po) reduce to three-dimensional inte- 
grals over the variables k, p, z , where z = cos 6, if 6 denotes the angle between the momenta 
k and p. 

At first we want to discuss the T = limit. An open question which occurs in the 
conventional NJL model is the choice of the cut-off for the second momentum integral 
in Eqs. (|23|) and (^51) over k. A very crude approximation presented in Ref. is the 
omission of the second integral by assuming that k = 0. In Refs. |13|JT5[1 the additional 
cut-off A was discussed. Ref. [|l3l assumes that A = Anjl and in Ref. [|15| upper and lower 
limits are determined from a calculation of in (9[1/Nc]. In the formulation we have 
chosen such a problem does not exist since the integrals are regularized in the separable 
approach by the proper treatment of the formfactors. The parameter Anjl in the cut-off 
formfactor ( ]26|) regularizes the integral over p. The upper limit of the k- integration is given 



by A = —pz + w4Anjl — — z'^) and runs between Anjl < A < 3Anjl- Note that in 
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solving (^3|) one has to check the integral limits for each term separately due to different 
combinations of formfactors partly hidden in the momentum dependent quark mass (^. 
Thus we have removed the ambiguity in regularizing the second momentum integration 
which occured in the previous approaches to the 1/Nc- expansion in the NJL model. 

The result for such a calculation in the T = limit is that for fixed model parameters 
the absolute value of the condensate is decreased by 20% compared to the Hartree- approx- 
imation. For comparison, a decrease of the quark mass due to the mesonic correlations in 
at T = has been obtained in Ref. [|l^]. This result can be understood qualitatively 
since the Hartree- contribution (p!0|) and the 1/Nc mesonic contribution (|23| ) to the quark 
condensate have opposite sign, the latter one being smaller in magnitude. 

Let us consider the finite temperature case. In order to compare the temperature be- 
haviour of the quark condensate in both models (Hartree approximation and Hartree ap- 
proximation with mesonic correlations) we have to fix the parameters such that the same 
values for the observables at T = are obtained, see Appendix A. The numerical evaluation 
of the final result for the quark condensate is shown in Fig. 5. Paying attention to the shape 
of the chiral phase transition, we observe that the inclusion of l/N^ mesonic correlations 
shifts the chiral symmetry restoration to lower temperatures when compared with the simple 
Hartree approximation. This finding is mainly due to the smaller qq coupling constant Kq 
for the model with mesonic correlations. 

VI. CONCLUSIONS 

In conclusion we have obtained the following new results: (i) a consistent l/N^.- expansion 
for the self-energy as well as for the propagator and a closed formula for the quark condensate 
up to the order l/Nc-, (ii) a finite temperature result for the l/N^ quark condensate within 
the Matsubara formalism, (iii) a consistent regularization of the two loop diagrams. 

The numerical evaluation for a NJL-type model shows that compared with the Hartree 
approximation mesonic correlations lead to a decrease of the absolute value of the quark 
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condensate at T = 0. After having compensated this effect of quantum fluctuations at T = 
by readjusting the model parameters (A, mo, Kq), the account for thermaly excited mesonic 
correlations shifts the onset of the chiral symmetry restoration to lower temperatures in 
comparison to the Hartree approximation. This behaviour seems to be closer to the recent 
results of lattice calculations where the condensate remains unchanged with temperature up 
to the chiral transition which occurs at T ^ 150 MeV for Nj = 2 p2 |. 



The inclusion of quark-antiquark correlations is of principal interest because the treat- 
ment of the medium in free quasiparticle approximation seems not to be appropriate at low 
temperatures. In contrast, in this region the mesonic degrees of freedom are expected to be 
relevant. This is supported by the fact that in the low temperature limit (T ~ 50 MeV) 
also other thermodynamic properties of a quark-meson system (e.g. the pressure) are dom- 
inated by mesonic contributions |T^. The presented l/N^.- expansion should be considered 
as a first step in including mesonic correlations. However, at temperatures where the chiral 
phase transition occurs, higher orders of the l/N^.- expansion may become important. 

Within the present approach, the treatment of the two-particle correlations was given 
in the usual pole approximation (^Dj) for the qq T-matrix. A next step in the evaluation of 
quark-antiquark correlations is the inclusion of the contribution of scattering states which 
will be considered in a forthcoming paper. In this way the account of the corrections due to 
two-particle correlations will be completed on the basis of the approach presented here. 
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APPENDIX A: TABLE OF INTEGRALS AND PARAMETER FIXING 



The polarization operators introduced in Eq. ([T9|) of the main text are defined as 



(Al) 



The temperature dependent meson masses are obtained from the solution of the Bethe- 
Salpeter equation 



1- J^(Po = M^(T),P = 0) = , 



(A2) 



where the polarization operators J'^{Pq) after evaluation of the Dirac trace and angular 
integrationare given by 



[1 - 2/(i?(q))] , 
[l-2/(E(q))] 



(A3) 
(A4) 



E{q) E\q)-{P,/2y 

The quark meson coupling constants introduced in Eq. (^) are evaluated in the rest frame 
of the qq pair (P = 0), where 



1 d 
2KoNf dP^ 



d 



Using Eqs. 



P2=MI 2KoNf 2PodPo 
and ([A4| ) , they are given by the following integrals 
1 /• , . o, . ^(q) 



(A5) 



J dq q^ g^iq) 
g-^^{M^) = ^Jdqq' g'{q) 



27r2 
1 



'2 1 



q 



l-2f{E{q)) 

l-2f{E{q)) 



(A6) 



The pion decay constant which we use for the parameter fixing at zero temperature is 
calculated by |llO| 

^cQTvqq f j^, ^,2 m^(q) 



fn 



dq q g{q) 



(A7) 



2^2 J ^ ^vh; E{q){E^{q)-M^/4) ' 

The model contains three parameters: the coupling constant Kq, the current quark mass rriQ 
and the range of the formfactor of the potential. We fix these 3 parameters to reproduce the 
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pion mass (M^ = 140 MeV) Eqs.(^ and the pion decay constant (/^ = 93 MeV) Eqs. 
( |X6| ) and ( |A7| ) and the quark condensate {— < qq >^^^ = 250 MeV) at zero temperature. 



The resulting parameter sets are given in Table 1. for both the Hartree approximation with 
and without the 1/Nc contribution from mesonic correlations, respectively. The parameters 
for the Hartree approximation are similar to those of the standard NJL- model, see Refs. 
ii|llJT||llj2l[ . 



APPENDIX B: MATSUBARA-SUMMATION 

The following formulae summarize the results of the one and two loop Matsubara-sums 
performed in this paper: 

dko I 1 - 2/(E(k)) 



2tt kl - E2(k) 2E(k) 

1 /i-2/(i?(p)) /(E(p))(i-/(£;(p))) 



(Bl) 



[dpo 1 _ 1 ( l-2f{E{ 

J 27r {pi - ^2(p))2 4^2(p) 2E(p) 



T 



+{i?(p) ^ -i?(p)} j , (B2) 
dko 1 1 + 2n(E^(k-p)) 



27r (fco - pof - El{\^ - p) 2E^(k - p) 

dko 1 
^{{ko-poY-El{],-mkl-E\].)) = 

/(E(k))+n(E^(k-p)) 



(B3) 



4(po + E^k - p) - E{k))E{k)E^{k - p) 
+{i5(k) ^ -Eik)}) + {E^ik - p) ^ -i?^(k - p)} , (B4) 



dpo 



27T {pi - E\p))ipo + E^ik - p) - Eik)) 
f{E{p)) - f{E{k)) 



2E{p){E^{k-p)-E(k) + E{p)) 



+ {^(P)--^(P)} • (B5) 
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Table caption 



Table 1. Sets of parameters for the Hartree approximation without and with the inclusion 
of mesonic correlations for a fixing scheme described in the text. 
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TABLES 



Ai")|")r()ximali())is 


A[M('V] 




/vo[MrV-2] 


Hartree 


660 


5.35 


9.45 


Hartreo+mcsonic correlations 


765 


4.4 


6.49 



TABLE I. 
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Figure captions 



Figure 1. The Dyson equation with the full self-energy. 

Figure 2. The Dyson equation in selfconsistent Hartree approximation. 

Figure 3. 1/A^c-expansion of the quark condensate. 

Figure 4. 1/Nc- approximation for the self-energy. The scalar- and pseudoscalar correla- 
tions are described by a RPA- type partial summation of bubble diagrams. 

Figure 5. The quark condensate as a function of the temperature in selfconsistent Hartree 
approximation (dotted line) and with inclusion of mesonic correlations (solid line). 
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